Bridge-mediated Donor- Acceptor Effective Coupling: Exact theoretical description. 
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An exact expression of the bridge-mediated donor-acceptor effective coupling, Hda, is derived. For 
systems represented by a tight-binding Hamiltonian with nearest-neighbor interactions, we show 
that |-ffda| 2 is equal to the product over all square of nearest-neighbor couplings divided by an 
appropriate product of level spacing of eigen energies of the Hamiltonian. Results of this calculation 
are compared to those obtained by perturbative approaches and some drawbacks of the latter are 
pointed out. 
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There are many physical, chemical and biological pro- 
cesses involving transfers of excitations between elec- 
trodes or from a donor (D) to an acceptor (A), through 
a mediated (environmental) composite bridge (B), that 
can be described formally by the reaction scheme, D* — 
B — A — ► D — B — A*. Such processes include, for in- 
stance, charge transfer, where charge are transfered from 
the donor to the acceptor through intermediate states 
connected them, or in energy transfer lily. LI . The quan- 
tity of interest for such donor-bridge-acceptor systems is 
the effective coupling H c i a between the donor and the ac- 
ceptor which carries the dominant distance dependence 
of charge transfer efficency, and that enters the transfer 
rate constant as a principal factor among others. Al- 
though the problem of computing the effective coupling 
is general to many processes, to be specific we will present 
our calculations within the framework of charge transfer. 
Recently, a number of experimental groups have reported 
measurements of the DNA charge transfer, and besides 
several theoretical models have been developed (see e.g., 
4] , 5] , more recent review papers 0,0, an d references 
herein). It is however disappointing that the results ob- 
tained so far using different approaches are not the same 
and even irreconcilable. It would therefore seem appro- 
priate at this time to put the theoretical treatment on a 
firm footing to analyze more accurately all essential in- 
gredients of the charge transfer phenomena, and among 
those the principal one - the effective D - A coupling Hd a 
is our concern in this paper. In what follows we derive 
an exact close analytical expression for Hda- It is not 
only of intellectual interest but also of relevance to gain 
further insight into the nature of DNA charge transfer, 
since comparing our exact expression and known in the 
literature results, one can estimate the accuracy and re- 
liability of various theoretical approximations. 

For a two-state system where the system is initially 
prepared in the donor state, the probability of finding 
the system in the acceptor state is given by, 



P d ^ a (t) = sin 2 



2h 



(1) 



ergies of the Hamiltonian describing the system. Us- 
ing the relation, lim sin 2 (xt) / x 2 = nt5(x), yields, 

t — >oo 

lim Prf_> a (t) ~ kd a t where hda is the rate constant given 

t — >OG 

by the Golden rule, k da = \ \H da \ 2 S{E + ~ EJ), where 
Hda = (E+ — E_)/2 is the effective donor- acceptor cou- 
pling of theproblem. For multi-state systems, several 

authors UlaHIHIISIHEim have used var ious 
approaches (perturbation theory, Lowdin's matrix parti- 
tioning technique, Green's function method) to derive an 
expression for the effective coupling. These approaches 
consist of mapping the multi-state eigenvalue problem 
onto the eigenvalue problem for the two-state system in 
order to determine the approximate splitting energies E± 
and H da as above. Accordingly, it is well established that 
the effective coupling between the donor and the acceptor 
can be expressed in terms of the Green's function of the 
bridge evaluated for the donor/acceptor energy [l2l llflj. 

To move further on smoothly and without loss of gen- 
erality, we consider a N + 2 system, {\n >} denoting 
the localized states on the the donor (n = 0), bridge (N 
sites with one state per site, denoted n = 1, 2, ••■ , N) 
and acceptor (n = N + 1). As pioneered by McConnell 
, we deal with the tight-binding (Huckel) Hamiltonian 
describing the charge transfer from a donor to an accep- 
tor through the bridge with site energies e„ and nearest- 
neighbor interactions v n between sites, 
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For this Hamiltonian (in which £0 = En+i = z), the split- 
ting energies E± can be obtained as roots of the equations 

s 



E±-£: 



i(E±) ± Vib(E±)} 2 + HEi 



(3) 



where, a(E) = v 2 G n {E) + v 2 N+1 G NN (E), b(E) = 



where E± (with E + > EJ) are the splitting eigen en- v 2 Gn(E) 



J N+ 



^NNiE), c{E) 



vlGn(E) 
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v 2 N+1 G NN (E), and G nm (E) =< n\(H - E) x |m > 
is the element of the Green's function restricted to 
the bridge described by its Hamiltonian Ho. Unfortu- 
nately, the solutions to these equations are too cumber- 
some. Thus, two main approximations are done in or- 
der to derive a useful expression of the effective cou- 
pling. First, for not too small bridge, it is tempt- 
ing to assume the following inequalities between the 
matrix elements of the Green's function, Gn(E±) ^> 
Gin(E±) and Gnn(E±) ^> G\n{E±). If these assump- 
tions are granted, neglecting the term c(E) in Eq.© 
leads to define an auxiliary energy Eq such that, Eq = 



e - v{G n (Eo) = e 



J N+1 



Gnn(Eo) (or equivalcntly, 



viGu(Eq) = vI !+1 Gnn{Eq))- Now, introducing this E 
back into Eq.© gives, E± w E ±viGin(E )vn +1 . Sec- 
ond, as it is still rather complicated to obtain Eq, an ad- 
ditional approximation consists in setting Eq ps e. Now, 
by the virtue of these two approximations, the effective 
coupling simply reads, 



\H. 



da I 



\viGi N (e)v NA 



(4) 



This is the well known formula, and the same result was 
arrived at by numerous methods 0, El , 0] . The su- 
perscript "(1)" refer to a sort of first order perturbation. 
It is so useful in practice that it is often used without 
care of whether the above associated approximations are 
still valid or not. For instance, Eq.© diverges at the 
resonance in the limit where the donor/acceptor energy 
coincides with eigen energies of the bridge. As shown by 
Zhdanov, in general Eq.© may overestimate the exact 
value by several orders of magnitude |l7j . Unfortunately, 
only a particular model for the electron potential energy 
has been investigated in , and this author did not pro- 
vide a general exact expression for the effective coupling. 

In this paper, we show that an exact expression for 
the effective coupling that does not suffer of above re- 
strictions can be worked out by using the long time limit 
of the time dependent solution of the Schrodinger equa- 
tion along with the Golden rule. To proceed, we solve the 
time dependent problem to determine the probability of 
charge transfer defined as, Pd-^ a (t) = | < 0| 1 3/(t) > | 2 , 
where \^f(t) > is the solution of the Schrodinger equa- 
tion, 

ih dm)> =HMt)> 

dt 

with the initial condition, |\l/(0) >= |0 >. If we denote by 
\4> n > and E n (with, n = 0,1, ■ ■ ■ , N + 1) the eigenstates 
and associated eigen energies of a given (TV + 2) x (N + 2) 
Hamiltonian, one can show that the transition probabil- 
ity reads as, 

N N+l 

P d ^a(t) = -J2 4<7v + i|n„|o> 

n— m=n-\-l 

(AE n:7n )t 



x <N+ l|n m |0> sin^ 
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where, AE n ^ m — E n — E m , is the eigen energy spacing 
and, Il n = 1 4> n ><(/>„ |, is the projection operator on 
the eigenstate \<f> n >. Equation J^J is a generalization 
of Eq.© to many state systems. As t — * oo, the lead- 
ing term in P<i^ a (t) will be associated to the smallest 
frequency \AE nm \/2?i corresponding to the smallest en- 
ergy spacing, IAEj+i^I, and defining hence the index I. 
Likewise, proceeding as done below Eq.©, we find for 
Eq.© that, lim P ( ;_ +a (i) ~ kd a t, with the rate constant, 

t — >oo 

kda = if \Hda\ 2 S{Ei+i — Ei), where the effective donor- 
acceptor coupling is given by, 



\H da \ 2 = - <iV+l|IL,|0> 

x <iv + i|n i+1 |o> (A£ M+1 ) 2 



(G) 



This expression provides a quite general and close for- 
mula for computing the effective coupling in various con- 
texts. 

Our main result can be formulated as follows. For 
a given Hamiltonian, compute the eigen energy spec- 
trum E n , sort the E n in ascending (or descending) order, 
and identify the index defined such that lAPz+i^l = 
min [|Ai? n!TO | ; n, m £ spectrum]. Next, determine the as- 
sociated eigenstates \4>i > and \(f>i+i >, and finally com- 
pute the effective coupling according to Eq.©. In ad- 
dition, an approximate expression of Eq.© can also be 
derived by using the perturbation theory or Green's func- 
tion technique for the projection operators H; and eigen 
energies E\. We expect that the |£/<j a | derived this way to 
be a better approximation because of the smallest eigen 
energy spacing A£y + i. 

Now, to illustrate the results in Eqs.© and ©, we 
derive explicitly the above expression for a tight-binding 
Hamiltonian as given in Eq. © . To this end, we consider 
the wave function of the form, \^f(t) >= Yln=0 c n{t)\n>, 
where c n (t) are the time dependent amplitude of the 
probability of charge at the nth site. This yields the 
following system of equations, 



ih 



r c + vici 



iti 



dc Q 
dt 

dc n 

ih = v n c n -i + e n c n + v n+1 c n+1 ; 1 < n < N (7) 

dt 

dcN+i 
dt 



VN+lCN + £Af + lCjV+i 



with the initial condition, c„(0) = 5 n Q. We define the 

Laplace transform f(s) = f£°dtf(t) e~ st of any function 
f(t). Laplace transforming Eq.© and solving the result- 
ing equation, we obtain the recurrence formula, 



IVN+l 



A„_ 1 c n+ i(a) + A„c n (s) = 11 (-^P) , (8) 



where A„ = (s + ie n /K)A n _i — (iv n /h) A„_ 2 with 
A_i = 1 and A = (s+ie /h). IiE n (n = 0, 1, ■ • • , N+l) 
denotes the eigen energies of H, we have: Aat + i = 



3 



JV+l 



Y[(s + iE n /h). As cn+2(s) = 0, wc obtain from Eq.©, 

n=0 



cjv+i(s) 



^JV+1 



N+l , . 

n 



JV+l 

n 

m— 1 



JV+l 



J] (s + iEJ„/7i) 



n=0 



After inverse Laplace transforming this expression we 
find that the probability Pd^a it) is given by Eq. (JSJ with 
the amplitude, 



<iv + i|n n |o>= 



JV+l 

n 

1=1 



N+l 

l[(En 

3=0 



(10) 



Now, using the formula in Eq.@) leads to the following 
expression for the effective coupling, 



JV+l 

n 



\H da \ z = 



JV+l 



(11) 



3=0 
jjiljil+1 

where the index "Z" is defined such that |Ai?; + i ;| is the 
smallest energy spacing of the eigen energy spectrum of 
the entire Hamiltonian. For this kind of Hamiltonian, 
the computation of the effective coupling is reduced to 
the determination of eigen energies of the system Hamil- 
tonian. 

As a direct application of the formula in Ea. (|ll() . we 
consider the following illustrative examples in which the 
energies of the donor and acceptor are both equal, £o — 
£ JV+l = £■ 

• System with 2 + 1 Levels: 

The simplest example which can worked out analyti- 
cally is the one where the bridge is reduced to a single 
site (i.e., N = 1). In this case, the Hamiltonian in Eq.J2J) 
reduces to a 3 x 3 matrix whose the eigen energies are 
given by, 



E = 
Ei = 
E 2 = 



s + ei- V(e-£i) 2 +4K+« 2 2 ) /2 



(12) 



e + e 1 + ^(e-e 1 y+4(vi+vi) /2 



It follows from this that the smallest energy spacing is 
AEi }0 with I = for E\ > e and AE 2 ,i with I = 1 for 
£\ < e. Applying the formula in Ea. ((TT|l . we get: 



\H da \ 



9 2 



1/ [AE , 2 A£i, 2 ] ; ei > e 
1/ [AE lt0 AE 2t0 ] ■ e x < e 



This is to compare with \H^\ = v\v 2 /\e\ — e\, obtained 
from using the approximate formula in Eq.Q). 
• Degenerate System with 2 + 2 Levels: 

An other interesting example where the formula in 
Eq.(@J fails to provide a finite value is, for instance, a 
bridge of size N — 2 with sites energies equal to the 
donor/acceptor energy, e± — e 2 = s. In this case, the 
eigen energies of the 4x4 Hamiltonian matrix are given 
by, 



E = e - A+/2 
E x = e- A_/2 
E 2 =e + X-/2 
E 3 =e + A+/2 



(14) 



where A± = J 2{v\ + v\ + u|) ± 1y/{v\ + v\ + vf) 2 - 4v 2 v 2 , 
with A+ > A_. Here, I = (or I = 2) corresponding 
to the smallest energy spacing AEifl (or AE 3y2 ). As a 
result, we find: 



\H da \ 2 = 



111 

■hv 2 v 3 



AEq >2 AE± t 2 AEo t 3 AE1.3 

Vivjv3 



4[(vi + v 3 ) 2 + > 



(15) 



• System with N + 2 Levels: 

As a last example, we consider the situation where the 
bridge of size TV is composed of sites with energies, e n = 
Eb and coupling, v n = v. In this case, the determinant 
A n for the eigen energies satisfies the recurrence relation, 



v„ - (e 6 - E)A n _ x + v 2 A n _ 2 = G ; n < N 
^jv+i - (e - E)A N + w 2 A/v-i = 



(16) 



with A_i = 1 and Ao = e—E. Introducing the variables, 
£(, — E = 2ux and a = (e^ — e)/w, one can show that, 



A„=t;" +1 [U n+1 (x)-2aU n (x)} 
Ajv+i = 2w(x — er)Ajv — w 2 A/v-i = 



n < N 



(17) 



where U n (cos8) = sin [(n + 1)0] / sin 8 is the trigono- 
metrical representation for the Chebyshev polynomials 
of the second kind that satisfy the recurrence relation, 
U n +i(x) - 2xU n {x) + £/■„_! (x) = 0. The characteris- 
tic equation for the eigen energies reads, Ajv+i = 0, i.e., 



U N+2 (x) - 2aU N+1 (x) + a 2 U N {x) = . 



(18) 



As a polynomial of degree N + 2 in x, this equation has 
N + 2 roots xj^V), fc = 1, 2, • • • , N + 2, which could 
be found numerically. To each root corresponds an eigen 
energy, Ek-i(cr) = £& — 2«x^(cr). One can convince our- 
selves when the eigen energies are arranged in ascending 
order (i.e., Eq < E\ < ••• < -Ejv+i, a choice for the 



(13) smallest energy spacing is E\ — E (i.e., x\ 



(JV) 



) so 



that 1 = 0. In this case, the effective coupling is given 



4 



by. 



\H A 



N+l 



n K> 



V 3 + l 



c j+l 



3=2 



As a check, Eg. (jT%)) reduces for N = 1 to, 8x 3 — 
8ax 2 + 2(cr 2 — 2)x + 2a = 0, whose solutions are: 



In contrast to Eg. i|19|) . this function diverges for a = 

and values a[ N) = 2 cos [kir/(N + 1)], k = 1,2, ■■■ , AT, 
(29) corresponding to the donor/acceptor energy equal to 
eigen energies of the isolated bridge, E^"^ 80 = Eb — 
2cos [A;7r/(Ar + 1)]. Numerical comparison of Ea. H22|) 
with the exact expression in Eq.lO shows that Eq.J23 
overestimates the value of the effective coupling. 



„(!) 



X [ L > = [\ a \ + Vo- 2 + 8j /4, 4 = \a\/2 and x\ 

yjcr 2 + 8 J /4. Using these roots in Eq.JTJJJ leads 
to Eq-CSJ). 

The evaluation of x^. (a) in explicit terms does not 
seem to be feasible. However, for a = or a = ±1, 
Ea. ((T%|l has the exact solutions (for k = 1, 2, • • • , N + 2), 



(0) 



A"7T 



(±1) = cos 



(k — l)?r 



7V + 2 



For comparison, the effective coupling obtained by using 
the approximate formula in Eq. Q is [13, lid] , 



Iff. 



2 N+1 v^a 2 ~ 4 



(a + V^ J3 4) JV+1 - (<r - V^4) 



JV+l 



To sum up, we have derived a quite general and close 
formula in Eq.JBJ for computing the effective donor- 
acceptor coupling for a given Hamiltonian. For a tight- 
binding Hamiltonian with nearest-neighbor interactions, 
Eq. © reduces to Eq. I|ll|) and the computation of the ef- 
fective coupling is reduced to the determination of eigen 
energies of the system Hamiltonian. The presented re- 
sults can be extended in several directions to the more 
general processes of electron transfer between two reser- 
voirs of states, exemplified by the molecular vibrational 
levels associated with the donor and acceptor sites con- 
nected by a molecular wire. It requires a set of different 
methods with specific approximations on certain scales of 
length and time but the coupling Hd a calculated in this 
paper is the indispensable entry point into such more re- 
(22) alistic study. 
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